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Abstract. We have studied various physical properties of magnetically deformed atoms and the associated 
matter, replacing the atoms by the deformed Wigner-Seitz (WS) cells at the crustal region of strongly 
\l ' magnetized neutron stars (magnetars). A relativistic version of Thomas-Fermi (TF) model in presence of 

strong magnetic field in cylindrical coordinates is used to study the properties of such matter. 
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Q^' 1 Introduction 

2 ■ From the observational evidence of a few strongly magnetized neutron stars, which are supposed to be the sources 
of anomalous X-rays and soft gamma rays, also called magnetars [T1I21I31I1] , the study of the effect of strong magnetic 
field on dense neutron star matter, including the crustal matter, both outer crust and inner crust regions of such 
'— compact stellar objects have gotten a new dimension. These exotic objects are also called anomalous X-ray pulsars 
^ . (AXP) and soft gamma repeaters (SGR). The outer crust of a typical neutron star in general, is mainly composed of 
■ dense crystalline metallic iron [^. The density of the upper edge of such metallic crystalline matter is ~ 10^ — 10 ''gm 
. cm~^, whereas the bottom edge is ~ lO^^gm cm~^, the matter consists of nuclei (also some highly neutron rich nuclei), 
surrounded by cylindrically deformed distribution of electron gas, makes the system electrically charge neutral and 
[ — drifted out free neutrons, if the density of matter is > neutron drift value. Since the density is high eneough, it is 
. therefore absolutely impossible to investigate the properties of such matter in material science laboratories, even at zero 
magnetic field. The observed surface magnetic field of the magnetars is ^ lO^^G, which is again too high to achieve in 
, the terrestrial laboratories. Also, it is quite possible that the interior field of such exotic objects can go up to ^ lO^^G 
(which can be shown theoretically by virial theorem) . If the magnetic field at the interior is really so high, then most 
of the physical and chemical properties of dense neutron matter should change significantly from the conventional 
neutron star (radio pulsar) scenario (see the recent article by one of the co-authors [5] for necessary references) . In 
[6] , the matter in the outer crust region of strongly magnetized neutron stars have been studied using Thomas- Fermi 
^ , approximation, in which the WS cells are assumed to be spherical in nature and arranged in a regular manner to 
H ' form dense crystal of metallic iron (~ BCC-type). Even though the magnetic field strength at the crustal region is 



slightly higher than lO^^G, it must change significantly most of the properties of dense matter, both in the outer crust 
and the inner crust regions of the magnetars [B][7]. It is believed that strong magnetic field can cause a structural 
deformation of the metallic atoms (see e.g., |S]) present in the inner crust region of a neutron star. The spherically 
symmetric structure of the atoms are destroyed and become cigar shape with the elongated axis along the direction 
of strong magnetic field. The atoms may even become almost an one dimensional string like object, i.e., needle shape, 
if the magnetic field strength is extremely high. In presence of ultra-strong magnetic field of strength > 4.4 x lO-'^^G, 
the application of TF model for spherically symmetric WS cells is not a valid approximation However, one can 
use TF model for sylindrically deformed WS cells because of ultra-strong magnetic field pi UnilTTlfT^ . In presence of 
ultra-strong magnetic field (^ 4.4 x lO^^G), the WS cells get magnetically deformed and become ellipsoidal in nature. 
In this article, for the sake of simplicity, we shall assume a cylindrical type deformation of the atoms in the inner 
crust region and use cylindrical coordinate system with azimuthal symmetry. In reality, to investigate the cigar like 
deformed atoms in presence of strong magnetic field one has to use prolate spheroidal coordinate system |13| . In future 
we shall present the problem related to structural deformation of atoms in a strong quantizing magnetic field using 
such coordinate system |13) . 
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In the present article we shall study the properties of inner crust matter composed of magnetically deformed 
metallic atoms. In section 2 we have developed the basic formalism and discuss the numerical results, whereas in the 
last section we have given the conclusions and the future prospect of this work. 



2 Basic formalism 

The width of the outer crust of a typical neutron star is ~ 0.3 km, the density of matter, which is assumed to be 
a dense crystalline structure of metallic iron is ~ 10^ — lO^^gm cm~'^, ranging from the upper edge to the bottom 
edge |14) respectively. In one of our previous work appeared iln [5], the properties of outer crust region has been 
studied. To investigate the properties of such dense exotic crystalline matter of metallic iron, we have replaced the 
outer crust matter by a regular array of spherically symmetric WS cells, with positively charged nuclei at the centre 
surrounded by a non-uniform electron gas. Whereas in the inner crust region, of width ^ 0.5 — 0.7km and density 
^ 10^^ — lO^^'gm cm~'^, we assume that since the magnetic field is high enough compared to outer crust region, the 
electron distribution around each nucleus (iron and also some neutron rich nuclei) gets deformed (the numerical values 
of widths and density ranges for various regions inside a neutron star strongly depends on the type of equation of 
state considered). They become cigar shape. In [7] we have studied the equation of states of inner crust matter with 
spherically symmetric electron distribution around each nucleus inside WS cells in presence of strong magnetic field. 
In this article for the sake of simplicity, we assume cylindrical type distribution of electron gas around each nucleus 
with the axis of each cylinder is along the direction of magnetic field and further assume azimuthal symmetry for 
all the cylindrically deformed WS cells. In this article, although we have considered magnetically deformed WS cells 
for electron distribution, the nuclei at the centre of these cells are assumed to be spherical in nature. We have also 
assumed that the magnetic field is not so strong to populate proton Landau levels and the magnetic (dipole) energy 
of neutron sector is negligibly small compared to the kinetic energy of these particles. 

To investigate various physical properties of deformed electron distribution in the inner crust matter, we start with 
the Poisson's equation, given by 

AttZc 

V^(t> = 47rene6l(r - r„) - — — 6l(r„ - r) (1) 

where Vn = 47rr^/3 = ^irr^A/i^ is the nuclear volume, rg = 1.12fm and A is the mass number of the nucleus [11] 
[T^ . In the cylindrical coordonate system, if one assumes the magnetic deformation of atomic nuclei in this region, the 
^-function associated with the contribution of protons within the nucleus has to be replaced by 0(r'„ — r)6{zn — ^), 
where in eqn.(l), Z is the atomic number of the nucleus, r„ is the nuclear radius, assumed to be spherical in shape, 
whereas in the above text, r„ and z„ given in the arguments of the ^-functions are respectively the radial and axial 
dimensions of the cylidrically deformed nucleus. However, in this article we have not assumed magnetic deformation 
of the nuclei. Here (/) is the electrostatic field, e is the electron charge and rie is the electron density, which because 
of assumed non- uniformity within the WS cells, is a function of positional coordinates (r, z). In our study, we shall 
consider the variation of <j)[r, z) within the cylindrical distribution of the electrons surrounding the positively charged 
nucleus. Therefore, in the Poisson's equation the nuclear part as shown in eqn.(l) will not contribute. Now in the 
cylindrical coordinate with circular symmetry, the above equation reduces to 

d^(j) 190 

It is well known that in presence of strong quantizing magnetic field, the number density of degenerate electron gas is 
given by 

where B is the constant external magnetic field, assumed to be acting along z-direction and is > Sc , where Be is the 
typical strength of magnetic field at and above which, in the relativistic scenario the Landau levels for the electrons 
are populated. For the sake of convenience, throughout this article we shall use natural units, i.e., h = c = 1. The 
critical strength is then given by Be — mj^/ | e | [Q, where me is the electron rest mass and | e | is the magnitude 
of electron charge. This critical strength may be obtained by equating the cyclotron quantum with the rest mass 
energy for electrons. We further assume that the matter is at zero temperature. In eqn.(3), pp is the electron Fermi 
momentum, v is the Landau quantum number, with Vmax, the upper limit of v. The upper limit will be finite at zero 
temperature and infinity for finite temperature. The factor (2 — S^q) takes care of singly degenerate v — Q state and 
doubly degenerate all other states with ^ 0. To study the properties of inner crust matter with deformed WS cells, 
we make Thomas- Fermi approximation, which is semi-classical version of Hartree approximation |15) . In this model, 
the well known Thomas-Fermi condition is given by 



{pp + ml + 2veBf/'^ ~ e(f> = constant (4) 
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where jie is the electron chemical potential, which is assumed to be constant throughout the WS cell. Hence one can 
express the Fermi momentum for electrons in the following form: 

PF = [{^Jie + e(|)f -ml-2ueBY/^ (5) 

Since the electrostatic potential (j) = <f)(r, z), the Fermi momentum pp for the electron is also a function of positional 
coordinates (r, z) within the cell. In principle one should use this exact expression for electron Fermi momentum in 
the equation for electron density (eqn.(3)) which in turn appearing on the right hand side of the cylindrical form of 
Poisson's equation (eqn.(2)). However, with this exact expression for pp, since the equation becomes non-linear in 
4)(r, z), it is absolutely impossible to proceed further analytically, even a single effective step. From the very beginning, 
therefore, one has to il technique to solve the Poisson's equation. Of course, with the numerical 

method within the limitation of the algorithm followed, more exact results can be obtained. However, in numerical 
computation of (t>{r, z), we only get a set of numbers, but the beauty of this model will be completely destroyed and a 
lot of interesting physics associated with the intermediate results of this problem will be totally lost. Therefore, to get 
an approximate analytical solution for (j){r,z), as a first approximation, we set the upper limit for Landau quantum 
number Vmax = and also neglect the rest mass of electrons, i.e., we put rUe = (since in the inner crust region 
the density of electron gas is high enough, the electron Fermi momentum will also be large compared to electron 
rest mass, therefore we expect that without appreciable error one can neglect electron mass in the above expression) 
in the expression for Fermi momentum pp (eqn.(5)). The approximation Vmax = is actually valid if the magnetic 
field is extremely high (~ 10^''G), perhaps is a valid approximation at the inner crust region for magnetars. However, 
to investigate some of the properties of dense electron gas within the cylindrically deformed WS cells to somewhat 
exact manner, later in this article, we shall use this approximate solution for 0(r, z) only, but do not use the first 
approximation Vmax = and rrie = to evaluate mathematical expressions for various physical quantities of the dense 
electron gas. However, later in this article we shall show that the upper limit t^rnax 

for the electron Landau quantum 

number is also a function of (r, z). Therefore to evaluate various physical quantities analytically in this region, first, 
one has to obtain an approximate solution for the Poisson's equation (with the values of I'max — and We = in 
the first approximation). To achieve our objectives, we use first the approximate form of electron Fermi momentum 
obtained from the assumption as mentioned above and its mathematical form is given by 

Pf » Me + (6) 

Next on substituting 

Aie + e(^(r,z) = V(r,2), (7) 
the cylindrical form of Poisson's equation reduces to 

^ + i^ + ^-A2V; (8) 

where = 2e^S/7r. Prom eqn.(8), it is quite obvious that under this approximation, the Poisson's equation reduces to 
a linear partial differential equation. To solve this equation analytically we use the method of separation of variables, 

given by 

ij{r,z) = R{r)L{z) (9) 
Substituting tp{r,z) from eqn.(9) in eqn.(8) and introducing a constant ^, we get 

(PR _^ 1 C^-R ^2^ _ g (•-|^Q^ 

f±-ie+x')L=o (11) 

where ^ is some real constant, independent of r and z but may change with the magnetic field strength and with the 
mass number and the atomic number of the type of elements present in this region. The solutions of eqns.(lO) and (11) 
are well known. For eqn.(lO), the solution is an ordinary Bessel function of order zero with the argument ^r, whereas 
for eqn.(ll), it is an exponentially decaying function of z. In the language of mathematics, the solution for ip{r, z) is 
then given by 

V'(r, z) = CJoi^r)exp \±{^^ + X'^^^^z'^ (12) 

where + and — signs are for z < or > respectively. We consider a convenient form of cylindrical coordinate system, 
such that z ~ plane is at the middle of the finite size cylinder. In that case we have to take positive sign for the 
upper half of the cylinder and negative sign for the lower half. Here C is a constant (again may change with the 
magnetic field strength and with the nuclear properties of the elements present in the inner crust region) and Jo{^r) 
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is the Bessel function of order zero. Now on the nuclear surface, at the centre of the WS cells, = Ze/r„ [TB], where 
Z is the atomic number and r„ = r^A^I'^ is the nuclear radius, tq = 1.12fm and A is the mass number. We consider 
eqn.(12) at various point on the nuclear surface. To evaluate the parameters C and ^ numerically, we do the following: 
put r = ar„ and z — /?r„ at a particular point on the nuclear surface and also consider r = /3rn and z — arn for 
another point, with + = 1. Then we have 



and 



From eqn.(14) we have 



V'(ar„,/3r„) = CJo(ar„C)exp [-{f + XY^^l3r,, 



C 



Ze" 



Jo(/3r„C)exp[-(e2 + A2)i/2ar„] 
Combining eqns.(13)-(14) we get 

Jo(ar„0 exp \-{e + X^Y'^Ptu] = MPr^O exp + X^^^arr, 



Hence we can write 



exp 



-ie+xY/\r.w~a) 



Jo(ar„^) 



(13) 
(14) 
(15) 

(16) 
(17) 



This is a highly transcendental equation for ^. However, it is possible to evaluate ^ numerically from this equation 
for a given magnetic field strength and for a given type of element, e.g., for metallic iron. To obtain ^ numerically we 
express eqn.(17) in the following convenient form 



e+x^ 



1 



In 



(18) 



To evaluate numerical values for ^, we now put (a) r and z values of the equator and poles of the nucleus, (b) 
r = r„/4, (c) r = r„/2 and (d) r — 3rn/4:, for four different positions on the nuclear surface and obtain four different 
sets of ^ as a function of magnetic field strength. For all these cases the ^-coordinates are obtained from the relation 
z = (r^j — r^)^/^. Because of the symmetry about z — plane, the choice of negative values for a and /? will give 
identical results. In fig.(l) we have plotted ^ (in MeV) as a function of magnetic field strength B, expressed in terms 

of critical magnetic field strength Bc'^\ for the specific values of r and z as indicated above by (a), (b), (c) and (d). 
For all these cases, the variations are insensitive for the low and moderate values of magnetic field strengths. This 
figure shows that beyond field strength 10^''G, ^ increases sharply for all the cases. It is also obvious from fig.(l) that 
the parameter ^ not only varies with the strength of magnetic field but also more strongly depends on the positional 
coordinates on the nuclear surface, giving one of the boundary conditionss. In fig. (2) we have plotted the same kind 
of variations for the parameter C with the magnetic field strength, expressed in the same unit as in fig.(l) and also 
considering the same positions on the nuclear surface as considered for fig.(l). For low and moderate magnetic field 
values it is also almost constant, then it falls abruptly beyond lO^^G and finally saturates to some constant values. 
Since beyond lO^^G, electrons within the cells occupy their zeroth or very low lying Landau levels, the quantum 
mechanical effect of magnetic field becomes extremely important in this region and as a consequence both ^ and C 
change significantly beyond this magnetic field value. Since the magnetic field is not ultra strong to distort atomic 
nuclei, we have therefore considered spherical nuclei of radii tqA^^^, at the z and r symmetric position inside the 
cylinder. From the solutions of eqns.(13) and (14), we have noticed that the functional form of the solution given by 
eqn.(12) does not change, but the unknown parameters C and ^ can have a large number of roots. As a result the final 
solution of Poisson's equation becomes degenerate (which has been shown in figs.(l) and (2)). As an hypothetical case, 
instead of spherically symmetric nuclei, we have replaced them by cylindrically deformed nuclei with axially symmetric 
nucleon distribution. The smaller cylinder (deformed nucleus) is coaxial with the bigger one and with identical z = 
plane. For the surface potential we have used eqn.(52) of Appendix A. In eqn.(52) we have replaced r^ax and Zmax by 
Vn and Zn, the maximum values of r and z for the deformed nuclei. Here ri and zi are the surface values of r and z 
coordinates at any arbitrary point on the nuclear surface. To make the nucleon distributions within the nuclei axially 
symmetric, we put 9i = 0. Of course, in such a geometrical configuration one of the geometrical parameter, e.g., either 
r„ or Zn has to be chosen arbitrarily, and the other one can be expressed in terms of the known one. Further, to choose 
one of the parameters, we assume that the nucleus is incompressible, as a consequence, the density will not change 
even if the geometrical configuration has changed. Then we have the simple relation r^z„ = SA/Arf,. In this expression, 
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if we choose one of the unknown arbitrarily (with the numerical value very close to the nuclear radius), the other can 
also be known. However, we have noticed that in this case the degree of degeneracy is even more than the spherically 
symmetric case. Along the curved surface for the nuclei, for ri = r„, we can have along the positive direction of z-axis, 
any value of z from to z„. Whereas on the plane faces, we have z = z„ and ri can have any value from to r„. All 
these points are on the cylindrically deformed nuclear surface and the potential on the surface is given by eqn.(52). 
Hence we may conclude that a non-degenerate solution can only be obtained, if and only if the electron distribution 
is spherically symmetric and the nucleus is also spherical with a common centre (concentric spheres) . Which we have 
studied previously [6 . However, we expect that if both the electron distribution and the shape of the nucleus at the 
centre are ellipsoidal in nature with the major axis (which is common for both the ellipsoids) along the magnetic field 
and two other axes (again common for both of them) are symmetric, the problem of degenerate solutions for electric 
potential inside WS cells will be removed. At present we are persuing this analysis. We believe that with this type 
of geometrical configuration, the electric field at the surface of the WS cell will also vanish, which is non-zero in the 
cylindrical case as discussed below.. 

Although the WS cells are overall charge neutral, at any point {rmax,z) on the curved surface and at any point 
if^Zmax) on the plane faces the potential ip{r^z) can not be a constant, here rmax and Zmax are transverse and half 
longitudinal dimensions. Which further means that the electric field at the surface (both radial and longitudinal 
components) can not be zero. This is a purely geometrical effect and such a deformed charge distribution exhibit 
quadrupole moment. The modified form of electro-static potential at any point on the curved face of the cylindrically 
deformed WS cell is given by 

il}{s,e) = CJo(C?-max)exp(-ylssin6l), (19) 
whereas on the plane faces, this potential is given by 

il){s,e) = CJo(Cscos6l)exp(-ylssin6') (20) 

where s = {r'^^x + z'^Y^'^ for the curved face and = (r^ -I- z^^^Y^^ for the plane faces, and A = (A^ + ^^-ji/a jg 
constant. Here the variable 9 is introduced to obtain the variation of z on the curved surface and also r on the plane 
faces. Therefore this 9 variable is not the conventional 9 coordinate used in cylindrical system. From eqns.(18) and 
(19), it is obvious that the potential can not be constant on the surfaces. The corresponding electric field on the curved 
surface is given by 

E = -^es + -^ee 21 
OS s o9 

= -2Cexp(-yl(s2 - rl^^^f'^){M^T^ax)Ae, + Ji(Cw)ee;) (22) 
Similarly the electric field at the plane faces is given by 

E = ^es + -^ee 23 
OS s o9 

^ -2C7exp(-ylz„„,)(Jo(e(s^ - zl,,Y'^)Ae, + .h{i{s^ - zl,,f'^)ier) (24) 

Which are obviously non vanishing on the surface of WS cells. However, neighboring cylindrical WS cells will interact 
electro-magnetically because of the non-zero values of electrostatic fields at the surfaces. Due to some kind of electro- 
magnetic induction there will be charge polarization on the surfaces of cylindrically deformed WS cells. As a result 
a number of charged (induced by the nearest neighbor WS cells) WS cells will form a bundle of charge neutral WS 
cylinders, instead of a single cylindrical cell 1 171. Again we expect that with ellipsoidal coordinate system, the non-zero 
electric field problem at the WS cell surface will also be solved. 
Now from the Thomas-Fermi condition, we have 

PF^[i^\r,z)~mlY'^ (25) 

where rrii, = (m^ -|- 2veBY/^ . With this exact expression for pp, the number density for electron gas can be expressed 
as 

n,(r, z) = 1^ (2 - 5,oW{r, z) - mlY'^ (26) 

i/=0 

This is obviously more exact than eqn.(3), where the value of Fermi momentum is taken from eqn.(6). Further, this 
expression shows that the electron density is a function of both r and z within the WS cells. Which actually justifies 
the assumption that the electron distribution inside each WS cell around the fixed nucleus is non-uniform. Now from 
the non- negative nature of p^p, we have 

_ %p'^{r,z) — nil ^ , \ 

^max — 71 7^ — ^raaxy ^ Zj j 
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The upper limit of Landau quantum number Vmax will therefore also depend on the positional coordinates r and z 
within and on the curved surface and the plane faces of the cylinder. In our model, the upper limit of Landau quantum 
number is zero at the deformed WS cell surface. Had the electron distribution at the crustal region been homogeneous, 
like electron gas in a highly conducting metal, with a background of positively charged nuclei at rest, the value of 
Vmax would have remain same at all points and for all the electrons depending on the density of electron gas and the 
strength of magnetic field. However, in our model, the electrons are distributed in an axially symmetric cylinder around 
a symmetrically placed positively charged nuclei of spherical in shape. The electrostatic potential changes with r and 
z accordingly following eqns.(9) and (12). From the nature of such variation, we found that Umax is maximum near 
the nuclear surface and vanishes at the WS surfaces (both curved surface and the plane faces). Therefore replacing 
both r and z by their maximum values Vmax and Zmax respectively, we have 

^maxiXmaxi ^max^ =0 (^^) 

Further, from the overall charge neutrality of the WS cells we can write 

Z = 2'Ke \ \ ne{r, z)rdrdz = constant (29) 

These two equations (eqns.(28) and (28)) are solved numerically to obtain Vmax and Zmax for various values of magnetic 

field strength for the inner crust matter with metallic iron only. Knowing both Vmax and Zmax , if we solve the equations 
Vmaxir.Zmax) = and Vmax{rmax,z) = numerically for r (ranges from to Vmax) and z (ranges from to Zmax), 
which are on the surface of the cylinder, one can generate the whole surface of the cylinder, including the plane faces. 

In fig. (3) we have plotted electron number density obtained from eqn.(26), in terms of normal nuclear density 
multiplied by 10'', as a function of radial distance from nuclear surface to the WS cell boundary for the magnetic field 
strengths lO"', 5 x 10^, 5 x 10^, 10** and 5 x 10'' times B'f^ ^ indicated by the curves a, 6, c, d and e respectively (In all 
the plots, we have used the values of C and ^ for the boundary condition on the equator and poles of the spherically 
symmetric nucleus. Although theoretically speaking the electrostatic potential inside the WS cells are degenerate with 
respect to the boundary condition on the nuclear surface, in practice, we have noticed from the numerical calculation 
that the variation with other points as the boundary is not so appreciable). These curves show that the electron 
number density is maximum near the nuclear surface and minimum near the WS cell boundary Tmax- This figure also 
shows that electron number density increases with the strength of magnetic field. In fig. (4) we have plotted the same 
quantity as in fig. (3) but against the axial distance z from the nuclear surface to the WS cell boundary (plane faces) 
indicated by Zmax- In this case also the variations are exactly same, qualitatively and quantitatively, as in fig. (3). 
The little qualitative difference is because of different types of functional dependencies. In fig. (5) we have plotted the 
upper limit of Landau quantum number Vmax as a function of radial coordinate from the nuclear surface to the WS 
cell boundary. In this figure the upper curve is for 10 x ^ the middle one is for 10^ x b'^^ and the lower one is for 
500 X Be ■ The value of Umax decreases with the increase in magnetic field strengths. It has been observed that beyond 
500 X Bc'^'' , Umax becomes identically zero throughout the WS cell. Further, the value of Umax for B < 500 x B^'^^ is 
largest near the nuclear surface and exactly zero near the cell boundary. In fig. (6) we have plotted the same kind of 
variation of Umax, but against the axial distance from nuclear surface to the cell boundary. The qualitative and the 
quantitative variations are exactly identical with fig. (5). These two figures show that the electrons are completely spin 
polarized in the direction opposite to the external magnetic field B at the cell boundary, including the plane faces, 

even for B < 500 x Bjf \ but beyond this value they are polarized at all the points within the cells. Although the 
variations along radial and axial directions are shown in these two figures, we expect that such polarized picture of 
electron gas exists throughout the cylindrically deformed WS cell surface, including the two plane faces. To explain the 
phenomenon of electron spin polarization within the WS cells, one has to solve Dirac equation in presence of strong 
magnetic field. It is then trivial to show that the Landau quantum number u is given by: 2u = 2n + 1 + a, where 

n = 0,l,2, is the principal quantum number and a = ±1 are the eigen- values of spin operator az corresponding to 

up and down spin states of electrons respectively. Hence it is quite obvious that all the Landau levels with u ^ are 
doubly degenerate (with two possible combinations of n and a), whereas zeroth order Landau level is singly degenerate 
with only possible combination is = is n = and a = —1. Which actually mean that in the zeroth Landau level 
spins of all the electrons are aligned in the direction opposite to the external magnetic field. This is a kind of electron 
spin polarization under the influence of ultra-strong magnetic field, occurring at the surface region of WS cells, and 
also throughout the cells beyond some magnetic field strength. This type of spin polarization will not be observed if 
electrons satisfy schrodinger equation, even if we introduce electron spin by hand. Therefore, the spin polarization is 
a purely relativistic effect. In fig. (7) we have shown the variation of Zmax with the strength of magnetic field B. This 
figure shows that the variation is almost insensitive for low and moderate magnetic field strengths but decreases almost 
abruptly beyond lO^^G, when most of the electrons occupy their zeroth Landau level, and finally tends to saturate 
to a constant value ~ lOfm. The variation of rmax with the strength of magnetic field is shown in fig. (8). The nature 
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of variation is more; or less same as that of Zmax- However, we have noticcid that for extremely large field strength, 
fmax — > 0, instead of saturation. This is a remarkable difference from its longitudinal counter part. It actually shows 
that in presence of extremely strong magnetic field the cylindrically deformed WS cells become more and more thin in 
the transverse direction. Wc therefore conclude that with the increase in magnetic field strength the radial contraction 
will be enormous compared to the axial one. From figs.(7)-(8) we have noticed that the variations are most significant 
beyond B = lO^^G. The reason is again because the electrons occupy only their zeroth Landau level in presence of such 
strong magnetic field, at which the qiiantiim mechanical effect of the magnetic field dominates. In these two figures 
we have taken Z = 26 and A = 56, which are tlie atomic number and mass number respectively for the deformed iron 
atoms. 

Next we calculate the different kinds of energies associated with the electron gas within WS cells. The cell averaged 
kinetic energy density of the electron gas is given by 

i/=0 
PF{r,z) 

dp,[{pl + miy/^-me] (30) 

where in the cylindrical coordinate system with azimuthal symmetry d^r = 2-jTrdrdz, with the limits r„ < r < rmax 
and r,i < z < Zmax and V = 7rr,^„,j,j.2zma,i., the vohimc of each cell. The pz integral is trivial, which gives an analytical 
expression for tlic local kinetic energy density i.e., at a particular point {r,z) within the WS cell. The factor 2 in 
the expression for V is for z-symmetry about z = plane. In fig. (9) we have plotted the variation of local kinetic 
energy density as a function of radial distance r from the nuclear surface to the cell boundary for various magnetic 
field values, keeping z = Q. This figure shows that the kinetic energy density increases with the increase in magnetic 
field strength. We have indicated the curves by a, b, c, d and e for the magnetic field strengths lOBjf \ 5 x 10'^Bjf\ 

5 X 10^Bi''\ 10^ si*^ and 5 x 10^ Bc*^ respectively. In fig.(lO) we have shown the same kind of variations along axial 
distance from the nuclear surface to one of the plane faces of the WS cell. The variation with magnetic field strength 
is again almost identical with fig. (9). Similar to the variation of electron number density within the cell (fig. (3) and 
fig. (4)), both from fig. (9) and fig. (10) it may be concluded that the local kinetic energy density for electron gas is 
maximum near the nuclear surface and minimum at the cell boundary. 

Similarly, the cell averaged electron-nucleus interaction energy per unit volume is given by 

Een = -yZe J d ry^—^^ (31) 

Analogous to the local kinetic energy density, here also one can obtain the local interaction energy per unit volume, 
Een{r, z) at a particular point within the cell. In fig. (11), the variation of the magnitude of electron-nucleus interaction 
energy per unit volume with the radial distance is plotted. In this figure wc have indicated the curves by a, b, c, d and 

e for the magnetic field strengths 10Bi''\ 5 x 10'^Bi''\ 5 x 10^Bi''\ 10"^ Bi''^ and 5 x lO^Bi""^ respectively. This figure 
shows that the magnitude of electron-nucleus local interaction energy increases with the increase in magnetic field 
strength. Which means that with the increase in magnetic field strength the electrons become more strongly bound 
by the Lorentz force of the form ev x B. In fig. (12) we have shown the same kind of variation along z-axis. Both the 
qualitative and the quantitative nature of variations with the strength of magnetic field are same as that of fig. (11). 

It is also obvious from figs.(9)-(12) that for a given magnetic field strength the kinetic energy density and the 
magnitude of electron- nucleus interaction energy per unit volume at a particular point, cither along axial direction 
or in the radial direction, inside WS cells are of the same order of magnitude. Wc do cixpect that the same type of 
variations will be obtained at all the points inside the WS cells. In fig.(13) we have shown the variation of cell averaged 
kinetic energy of electron gas (solid curve) and the magnitude of cell averaged electron-nucleus interaction energy 
(dashed curve) with the magnetic field strength. This figure shows that both the quantities are increasing function 
of magnetic field strength. However, for low and moderate field values, the variations are not very much sensitive, 
particularly for the kinetic energy part of electron gas. 

Next we consider the cell averaged electron-electron direct interaction energy density, given by 

Etr = J d'rn^ir, z) j d\'n^{r\ _ _ (32) 

Assuming r as the reference axis, we have (r — r')^ = + r'^ — 2rr'cos6', where Q is the angle between r and r', 
assumed to be on the same plane. Then d^r' = r' dr' dOdz' , with < < 27r. In this case the z integral has to be 
broken into two parts, one with limit —Zmax ^ z <'z and the other one with the limit 'z < z < +Zmax- The value of 
z is not easy to evaluate in the region between z-axis and r-axis. With 9 symmetry, for the sake of simplicity we put 
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z = r„ and expect that the error will be nominal. Now to obtain electron-electron direct interaction energy one has 
to evaluate the five dimensional integral as shown in eqn.(32). None of them can be obtained analytically, hence it is 
necessary to follow some numerical methods. Even the 6 integral can not be obtained analytically. One can express 
the 6 integral in the form of an elliptical integral of first kind, given by 

lecir, r\ z, z') = jj d^^^—^A-^ (33) 

where K = 4rr'/[(r + r')^ + (-^ ~ ^')'^]- The direct part is then given by 

Eii^ = ^e'J d'rn,{r,z) J r'dr'dzn,{r' , z') ^^^^ _ J,^(r, (34) 

where the factor 4 is coming from the angular integral over 6 from to 27r. Let us now consider the elliptic integral 
(eqn.(33)) on z = 0-planc. In this case both z and z' are zero and the factor K = 1. Then it can be shown very easily 
that the integral given by cqn.(33) will diverge at the lower limit. To avoid this unphysical infinity wo put a lower 
cut-ofF (infrared cut-off) S, which will now the lower limit for 0. The physical meaning of non-zero lower limit for the 
0-integral is that the two electrons under consideration can not be at zero distance from each other on the arc of a circle 
whose centre is same as that of the nucleus. The infrared cut-ofF S which is a measure of minimum angular distance 
between two neighboring electrons must necessarily dcpc^nds on the minimum possible linear distance between them 
and also on the average radial distance from the centre of the nucleus. From a very elementary geometrical construction 
it can be shown that 

5=- (35) 
r 

where s is the arc length, or the distance between two neighboring electrons on the circular arc. Since s is infinitesimal 
in nature, we can approximate it by a straight line of length s, which is the length of the cord connecting two points 
occupied by two neighboring electrons. Since s is the minimum possible linear distance between two electrons, we can 
express it in terms of electron density near those points, given by 

s ~ n-V3 (36) 

The coll averaged electron-electron direct energy has been obtained by evaluating the multi-dimensional integrals 
numerically. 

The electron-electron exchange energy corresponding to the ith electron in the cell is given by 

EiT' = -yE/ d'rd'^' [(^ _ +\z - ,,)2^i/2 Mr,z)Mr',z')Mr,z)Mr',z') (37) 

where ^i(r, z) is the spinor wave function of Dirac equation in cylindrical coordinate in presence of strong quantizing 
magnetic field, and ^p{r,z) = il)^{r,z)'yQ, the adjoint of the spinor and 70 is the zeroth part of the Dirac gamma 
matrices 7^ in cylindrical coordinate system. We have evaluated the cell averaged exchange energy using Dirac spinors 
in cylindrical coordinate. An elaborate discussion is given in the Appendix. 

The kinetic pressure of non-uniform electron gas within the WS cell is given by 



The Pz integral is very easy to evaluate analytically and is given by 

47r2' 



(38) 



;y=0 



_„>{?£±fe4±I=^}] ,39) 

which gives the local pressure at (r, z). This equation shows that the electron kinetic pressure also changes from point 
to point within the WS cells. In fig. (14) we have shown the variation of kinetic pressure for the non-uniform electron 

gas with the radial distance within the cell. Curves a and h are for B = 10 X bI:^ and B = 500 x 5,^"^ whereas upper 

curve and lower curve as indicated by c (almost identical) are for B = 5000 X b'^'^ and B = 10000 x B'f'^ respectively. 
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and 

(e) 



In fig. (15) the same kind of variations are shown along z-axis. In this figure the curves for B — 5000 x B, 

B = 10000 X Be are almost identical and indicated by single thick curve c. For B = 10 X B^''^ and B = 500 x B, 
the curves are indicated by a and 6 respectively. These two figures show that the kinetic pressure is maximum near 
the nuclear surface and zero at the cell boundary. The variation with magnetic field strength shows that the non- 
uniform electron gas becomes softer for higher magnetic field within the cells. In fig. (16) we have shown the variation 
of cell averaged electron-electron direct interaction energy (solid curve) and the corresponding kinetic pressure (dashed 
curve) with the magnetic field strength. Both the quantities are monotonically increasing function of magnetic field 
strength. Since in TOV equation for neutron stars, in the equation of state the cell averaged electron gas kinetic 
pressure is used, hence we can conclude that since electron gas becomes harder in presence of strong magnetic field 
in this particular region, the width of inner crust region will increase with the increase in magnetic field strength. 
In fig. (17) we have shown the variation of cell averaged electron-electron exchange interaction energy, assuming the 
extreme case when electrons occupy only their zeroth Landau level (dashed curve) and the most general one, when 
electrons can have all possible Landau levels (solid curve). The figure shows that for ^ 0, the exchange energy is 
oscillatory for low and moderate values of magnetic field strength. This oscillatory phenomenon is identical with the 
observed De Haas- van Alphen oscillation observed in Landau diamagnetism. The exchange energy becomes extremely 
small when the electron Fermi momentum suddenly becomes zero for some value of magnetic field strength and then 
rises sharply upto a certain magnitude of magnetic field strength and finally decreases to zero, when the magnetic 
energy dominates over the matter part. In the case of = 0, the nature of variation is almost identical, except the 
oscillatory nature at low and moderate magnetic field region. 

For the sake of completeness, we have obtained the classical form of electron-electron Coulomb potential energy 
and the corresponding electron-nucleus interaction energy within a WS cell. The electron-electron interaction energy 
is given by (an elaborate discussion is given in the Appendix) 



Epf, = 



cJi{kr 



[k{2 



^5 ['max'- 

n) - (1 - exp{-fc(z„ 



-rn)})] 



and the corresponding electron-nucleus interaction part is given by 



' max^rriax 



(40) 



(41) 



where 



/ = 



fz^ _L^2 Nl/2 , 2 log f (^'^^ax + rliax)^ ^ + 

K'^max ^ ' max) ^rnax ^ ' max -'-"-'6 I 

\ 'max 

(r^+r^ ,^2 log ( (^n+rmaxV^^ + ^« \ 

Vn<'max) ' n ^ ' max \ I 

\ 'max J 

iz^ -^r^S^I'^Z I log i ^ ^nY'"^ + Zmax \ 

y/^max ' ' n) -^max ^ ' ,i ^'^6 I ^ ) 



+ ^[2^/M + r^log(2^/^ + l)] 

(42) 

The detail derivations for electron-nucleus Coulomb interaction energy is also given in the Appendix. In fig. (18) we 
have shown the variation of electron-electron Coulomb energy density (solid curve) and the corresponding magnitude 
of electron-nucleus classical interaction energy density (dashed curve) with the strength of magnetic field. Both the 
curves show that the classical interaction energies are also insensitive to the strength of magnetic field in the low and 
moderate regions, and both of them are affected significantly beyond lO^^G. The reason is again that the electrons 
occupy only their zero-th Landau level. However, unlike the quantum mechanical cases, here at high magnetic field 
region both of them become extremely small. Moreover the overall magnitude does not change by an order of magnitude 
within the range of magnetic field considered. Finally to show that the cylindrically deformed WS cell structure of 
the inner crust matter of strongly magnetized neutron star (magnetar) is energetically favorable over the spherical 
structure, in fig. (19) we have compared the energy per electron for these two possible type of WS cell structures for 
various values of magnetic field strength. In this figure solid curve indicated by sph is for spherical cell structure and 
the dashed one indicated by cyl is for the cylindrically deformed WS cells. The energy per electron plotted along left 
y-axis is for the spherical case whereas for the cylindrical case the same quantity is plotted along right y-axis. This 
figure shows that for cylindrical case the total energy per electron is about one order of magnitude less than that 
of spherical case. Further beyond B ^ lO^^G, energy per electron for the cylindrical case becomes several orders of 
magnitude less than the spherical case (see also [T5]). 
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3 Conclusions 

In this article we have investigated various physical properties of non-uniform electron gas assuming cylindrically 

deformed atoms of metals, in particular the metallic iron at the inner crust of a strongly magnetized neutron star. 
Because of extremely strong surface magnetic field of magnetars, we have assumed a cylindrical type deformation of 
the atoms, which are subsequently replaced by WS cells with the same kind of geometrical structure. The longitudinal 
axis of all the cylinders are along the direction of magnetic lines of forces. The curved surfaces of these cylinders are 
therefore parallel to the boundary surface of the neutron stars in the region far away from the magnetic poles, where 
the magnetic lines of forces emerge almost perpendicularly with the surface (polar cap). We have studied various 
physical quantities for the dense electron gas in the inner crust region. Wc have investigated the variations with 
magnetic field strength for cell averaged quantities and also the spatial variations for constant magnetic field. We 
have noticed that the transverse dimension of a cylindrically deformed WS cell becomes extremely thin in presence of 
ultra-strong magnetic field. In this work we have also compared the total energy per electron with the spherical case 
and found that the cylindrically deformed WS structure in the inner crust region is energetically more favorable over 
the spherical case in presence of ultra-strong magnetic field. Although we have considered cylindrical type deformation 
for the iron atoms in this region, it is expected that the atoms become cigar shape in presence of strong magnetic field. 
The present investigation, assuming cylindrically deformed atoms is therefore an approximate model calculation. In 
our future study, the properties of neutron star inner crust matter with cigar shape atoms in the metallic crystal in 
presence of strong magnetic field will be investigated. 



4 Appendix 

A. Coulomb energy of electron gas: 

1. Electron- nucleus Coulomb interaction energy: 



iJen = -Ze I ^ (43) 



Where the integral is over the whole volume and s = (r^ + z'^Y^^. The charge element within an elementary volume 
dv is given by 

, ^ dv rdrdz 

dq = Ze- = Ze (44) 

" ^' max'^max 

Hence 

^2g2 ^ 
^en = -^^^ / rdr / ( 2 + ^2)1/2 (45) 

The double integral can very easily be evaluated and finally we get eqn.(40) as given in the text. 
2. Electron-electron Coulomb interaction energy: 

Where s = [{zi — Z2Y + (ri — r^YY^"^- In this case the charge elements dqi and dq2 are at {ri,9i,Zi) and (r2,02,^2) 
respectively. The factor 1/2 is to take care of double counting. Prom a very simple geometrical construction it can 
very easily be shown that 

s = [(zi - z^f + (r\ +rl- 2rir2 cos(^^i - ^^2))] ^'^ (47) 
Substituting for dqi and dq2, with the definition as given in eqn.(42), we have 

Eee = 7: . o A o — / ^idri / r2dr2 / dzi / dz2 / d9i / d62 



2 47r2^4 ^2 



1 



[(^1 - Z2Y + (r-f +4- 2nr2 cos(^i - $2))^^ 



(48) 



To evaluate the integrals, we use the identity 
1 



« [(^1 - Z2r + {rl +rl- 2nr2 cos(^i - ^2))] 



1/2 



- exp{im{9i - 62)} Jmikri)Jm{kr2)exp{-k{z> - z<)}dk (49) 
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To obtain electron-electron Coulomb energy, we first evaluate the potential at {ri,6i,zi) due to a charge element 
dq{r2,92, Z2). This is given by 

(l){ri,0i,zi) = / (50) 

Using the identity as defined above (eqn.(48)) and integrating over Z2, within the range r„ < 02 < Zmax, we get 



{ri,Oi,zi) 



V] / -j- r2dr2 I d92Jm{kri)Jm{kr2) 

- Jo « Jrn Jo 



(2 - exp{-k{zi - r„)} - ex.p{-k{zmax - zi)}) exp{im{Oi - 62)} (51) 



Integral over 62 gives 27r^mo, which means only m = term of the series in the above identity (eqn.(48)). Hence, we 
have 

Ze f°° dk 

</>(ri,6'i,2;i) = ^ / — / r2dr2Jo{kri)Joikr2) 

'^max^max JQ Jr„ 

(2 - exp{-fc(^i - r„)} - ex.p{-k{zmax - zi)}) (52) 
Further, to evaluate the r2 integral within the range r„ < r2 < rmax, we use the relation 

^ xr{^) ^ {\m\+2i + i)r{f + i) ^ 

rdrJmir) = pf^^^^) ^ p(m+4 , J|„|+2;+i(a;) (53) 



f 

Jo 



Since m = 0, we have in the denominator of the above expression r{m/2) = 00. Whereas the mimcrator of the first 
term of I series in the above relation (i.e., for / = term) also contains a r{m/2) = 00 term. These two terms will 
cancel each other and we get non-zero contribution for Z = only with m = 0. To evaluate r2 integral we decompose 
it in the following form as given below 

nax P'^max P'^n 

dr2 = / dr2 - / dr2 (54) 



Then using eqn.(52) with I = and m = 0, it is possible to evaluate r2 integral analytically. Then we have the Coulomb 
potential due to the charge element dq{r2,02, Z2) 



</>(ri,^l) = -2-— / t| [rmaxJl{krmax) - rnJl{krn) 

'^maxZmax Jo 



' max^max 

Jo(fcri) [2 - exp{-fc(^ - r„)} - exp{-k{zmax - zi)}] (55) 

Hence it is also possible to evaluate the components of electric field, — —dip/dzi along z-dircction and £",. = 
—d(f)/dr\ along the radial direction. It is trivial to show that these components are non-zero on the WS cell surface, 
including the plane faces. The Coulomb energy is then given by 

Eee = -^-z — 2 / (li{ri,zi,ei)dq{ri,zi,ei) (56) 

^ ^^'^maxZmax J 

The 9i integral will give 2tt, the zi integral is very easy to evaluate and for ri integral wc use the Bessel integral 
formula as given in cqn.(52). Then we have the electron-electron Coulomb interaction energy 



■ J \j'maxJl(kl'max) ~ ^n>^l(^^n)] [^(•2max ~ fn) ~ (1 ~ 6xp{ — fc(^max ~ ^n)})] (5'^) 

B. Electron-electron exchange energy: 

In the cylindrical coordinate system with the external magnetic field along z-axis, which is also the symmetry axis 
of the cylinder and considering the gauge A(r) = B{—y/2, x/2, 0), we have the Dirac equation corresponding to upper 
component 4>\{r): 

E^-m' + 2Xk+-^ + -^ + ^- epA Mr) = (58) 

op'' p op oz^ J 

where k = eB/2 and az4>{f) = \4> with A = ±1, eigen-values for spin-up and spin-down states respectively. Defining 
j3\ = E"^ — w? -\- 2Xk, the above equation can be written in the following form. 



15 5' ,2 2 

op^ pop OZ'^ 



4>x{p,z)=Q (59) 
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writing the solution for the upper component in the separable form: (j)\{p,z) = fx{p)exp{ip;^z), we have from the 
above equation 

\o 1 ^ ,22] 

op'' p op 

where I3x is replaced by P\— p^. The solutions are given by 



fxip) = 



(60) 



fx{p) = exp 



9x{t) 



(61) 



where t = rj^ = kp^. Substituting f\{p), we have 

9x = exp{ip^z) exp ( -- j L„(t) 
where L is the linear dimension along z-axis, Ln{t) is the well known Laguerre polynomial, 



(62) 



N 1= 



1/2 



(63) 



is the normalization constant. Considering the effect of spin into account, we have 



gx=+i = Ln [ ^eBp' 



(64) 



and 



9x=-i = Ln Tif^Bp- 



(65) 



are the up-spin and down-spin states respectively for the upper component, with the corresponding energy eigen-values 
E+ = {pI + + 2neBf'^ and E+ = {pi + w? + 2{n + VjeBf/^. 
The spin-up and spin-down spinor states are then given by 



k /_ m\ 



exp(ip^2;)exp(-|) 



Ln \ 



P;,Ln/iE + m) 
-2ik^/H^/^/{E + m)L'„/ 



(66) 



and 



il>^ip,z) = 



k ( m\ 
2Ei'+E) 



1/2 



exp(ipz2:)exp(-|) 





-2ik^/H'/y{E + m){L'„-Ln) 
-p^Ln/{E + m) 



Substituting for up-spin state, we have 

l^]{M{r') = ^ (1 + I) cxp{-zp,(z - z') - \{t + t')} 

Pi \ L (t)L (t') - ^^^L 



2ttL 
1 - 



(67) 



(68) 



Similarly, 



2ttL 
1 



(1 + ^) exp l^-ip,{z - z') - ^(i + i')| 



{E + m)2 



Akitt')^/'^ 



(69) 



Identical expressions can also be obtained for jth type particles, on which there is a sum. To obtain the exchange 
energy, we replace the sum over j by the integral over momentum p^ of jth particle, i.e., Ylij ~^ J dp'^. Assuming r 
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as the reference axis, we can write J (PxcPx' ... — 27r J pdpdzp'dp'dz'dd' .... Replacing the variables zf and p'^ by two 
new variables, z and Pz, where z — z' — z and = pz — p'^, we have the z integral 

is the modified Bessel function of order zero. The z-integral simply gives L. The integral over P^ can be obtained 
analytically using the following integrals (eqns.(71)-(72)): 



where X ^\ p — p' \ and 



2pF 



1 



/ = / KoiPz I X \)dP, = —— / Koiw)dw (72) 
Jo \ ^ \ Jo 

where a — 2pp \ X \ and pp is the electron Fermi momentum. The value of the integral on the right hand side of 
the above equation is obtained from standard mathematical hand book [T5] and is given by 



/ Ko{t)dt = -{j + \ogx/2) xY^jr 
Jo 



{x/2) 



2k 



(fc!)2(2fc + l) 
(a;/2)2fe ~ (x/2)2'= 



(fc!)2(2fc + 1)2 (fc!)2(2fc + 1) " 

(l + l/2+... + l//c) = x/(x) (73) 

where 7 (Euler's constant) — 0.5772156649. The rest three integrals, over p, p' and 9 are evaluated numerically. The 
Laguerre polynomials Ln{x) and their derivatives L'^{x) are obtained numerically from the recursion relations and 
derivative formulas from the reference as cited above and finally obtained the exchange energy as a function of magnetic 
field strength as defined in eqn.(36). 
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Fig. 1. The variation of separation variable ^ with the magnetic field strength. Curve (a): r and z are on the equator and poles 
respectively, curve (b): r = r-n/^-, curve (c): r = rn/2 and curve (d): r = Zrn/^- For the curves (b),(c) and (d), z = (r^ — r^Y^'^ . 
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Fig. 4. The variation of electron number density expressed in terms of normal nuclear density with radial distance z in Fermi. 
The curves indicated by a, 6, c, d and e are for S = 10, 5 x 10^, 5 x 10^, 1 x 10* and 5 x 10* times Bc'^' respectively. 
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Fig. 5. The variation of upper limit of Landau quantum number with the radial distance in Fermi. Upper curve is for B = 
10 X Bi''\ middle one is for 10^ x si"^ and the lower one is for 5 x lO'^Bi'^K 
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Fig. 6. The variation of upper hmit of Landau quantum number with the axial coordinate. Upper curve is for S = 10 x Bi''\ 
middle one is for 10^ x Bo*^' and the lower one is for 5 X 10^ Sr'. 
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Fig. 7. The variation of Zmax with the magnetic field strength. 




Fig. 8. The variation of rmax with the magnetic field strength. 
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Fig. 10. The variation of kinetic energy density with the axial coordinate in Fermi. The curves indicated by a, b, c, d and e 

are for B = 10, 5 x 10^ 5 x 10^, 1 x lO" and 5 x lO" times B^"' respectively. 



24 Please give a shorter version with: \authorruiming and \titleruiming prior to \maketitle 




4.3 4.4 4.5 4.6 4.7 

r (f m) 



Fig. 11. The variation of the magnitude of electron nucleus interaction energy density with the radial distance in Fermi. The 
curves indicated by a, b, c, d and e are for B = 10, 5 X 10^, 5 X 10^, 1 X 10'' and 5 x 10* times Bi"^ respectively. 
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Fig. 12. The variation of the magnitude of electron nucleus interaction energy density with the axial coordinate in Fermi. The 
curves indicated by a, b, c, d and e are for B = 10, 5 x 10^, 5 x 10^, 1 x 10* and 5 x 10* times Sc*^' respectively. 
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Fig. 14. The variation of kinetic pressure of non-uniform electron gas within WS cell with the radial distance in Fermi. Curves 
a and b are for B = 10 x Bi"^ and 5 x 10^ Bc*^' respectively, whereas upper and the lower curves indicated by c are for 5 x 10^ B^*^' 
and 10* X Bi^^ respectively. 
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Fig. 15. The variation of kinetic pressure of non-uniform electron gas within WS cell with the axial coordinate in Fermi. Curves 
a and b arc for _B = 10 x si"'' and 5 x IQ^B'f^ respectively, whereas upper and the lower curves indicated by c are for 5 x IQ^B^"^ 
and lO** X Bi^^ respectively. 
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Fig. 16. The variation of coll averaged electron-electron direct interaction energy in ev (solid curve) and the cell averaged 
electronn kinetic pressure (dashed curve) with the strength of magnetic field. 
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Fig. 18. The variation of cell averaged electron-electron Coulomb energy density (solid curve) and electron-nucleus Coulomb 
energy density (dashed curve) with the strength of magnetic field. 
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Fig. 19. Variation of total energy per electron in MeV with the strength of magnetic field. The solid curve is for the spherical 
case with energy per electron plotted along left y-axis and the dashed one is for the cylindrically deformed case with energy per 
electron plotted along right y-axis. 



